The Casimir energy between nanostructured gratings of arbitrary periodic profile 
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We study the lateral dependence of the Casimir energy for different corrugated gratings of arbi- 
trary periodic profile. To this end we model the profiles as stacks of horizontal rectangular slices 
following the profiles' shape and evaluate numerically the Casimir energy between them for different 
relative lateral displacements of the two corrugated plates. We compare our results with predictions 
obtained within the proximity force approximation (PFA). At comparable separation of the corru- 
gated plates and geometric parameters, we find a strong dependence of the Casimir energy on the 
shape of the corrugation profiles. 
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I. INTRODUCTION 

Two surfaces in quantum vacuum attract each other 
due to the fluctuations of the electromagnetic field pre- 
vailing in the vacuum state as predicted by H. Casimir in 
1948 [1 . The Casimir force scales with the inverse of the 
fourth power of the surfaces' separation and can take on 
significant values between objects in the submicrometer 
range. Indeed, the dynamics of the electromagnetic field, 
including the vacuum field, strongly depends on its exter- 
nal boundary conditions, and the variation of its spectral 
density when the field is confined leads to differences in 
the vacuum radiation pressure resulting in a net force. 

While the Casimir force is negligible at macroscopic 
scales it can become comparable to electric forces at 
submicron distances. In microelectromechanical systems 
(MEMS) it may produce sticking and adhesion of the sys- 
tem's components but it may also be used to actuate 
the system [31 14] . Modern measurement techniques have 
allowed for accurate determination and measurements of 
the Casimir force for different configurations and geome- 
tries. Casimir's original configuration of two fiat parallel 
plates has been experimentally realized in [5]. Most of 
the measurements have however been performed between 
fiat and spherical mirrors [H [6UTT] , and some have em- 
ployed cylinders [121 US] • 

Deviations from the two plates geometry are often ac- 
counted for in calculations by using the proximity force 
approximation (PFA) which amounts to adding up local 
contributions to the Casimir force at different distances 
as if those were independent of each other [14] . Exact cal- 
culations have shown considerable deviations from PFA 
for the plane-sphere geometry [15] |16] and when nanos- 
tructured plates are used jEl HI] • In particular the nor- 
mal Casimir force is affected by the gratings geometry, 
and the change cannot be reliably estimated within PFA 
as soon as the different length scales of the problem (cor- 
rugation period depth a and relative distance L) are 
of the same order. This is shown by the experiment re- 
ported in [T9^ for the interaction of a Au sphere with a 
Si grating with deep rectangular trenches. Exact calcu- 
lations within a path integral approach of this geometry 
for perfectly reflecting plates [20^ confirm the invalidity of 



PFA but lead to too large a prediction for the force. Tak- 
ing into account in the exact calculations the interplay 
between the surface geometry and the optical properties 
of the materials leads to good agreement between the ex- 
periment and theoretical predictions [17], as confirmed 
later on by another experiment using a shallow-trenched 
Si surface and comparing with independent theoretical 
calculations [18 . More recently these methods have been 
complemented by another approach based on a decom- 
position of the field onto the modes of the structures [2j . 
See also a recent review on Casimir forces in structured 
geometries by Rodriguez et al [22] , 

The use of nanostructured plates instead of flat sur- 
faces also allows for the existence of lateral Casimir forces 
acting tangential to the surface and of Casimir torques. 
They arise when the translational invariance is broken, 
either due to the anisotropy of the material [23] or to 
surface structures [24]. Lateral Casimir forces were first 
calculated for scalar fields and perfect boundaries im- 
printed with one-dimensional (ID) sinusoidal corruga- 
tions [25] [26], where by ID we mean corrugations fol- 
lowing a single transverse direction. The case of elec- 
tromagnetic fields and non-dissipative metallic surfaces 
was studied within the scattering approach in [27] and 
later on applied to Casimir-Polder interactions between 
a metallic nanostructure and Rb atoms [28]. The Casimir 
interaction involving single objects of elliptical shapes 
have been studied theoretically in [29H3T] and experimen- 
tally in [13 . 

Experimental evidence of the lateral Casimir force was 
given for sinusoidal corrugations imprinted on a sphere 
and a plate, with corrugation periods d much smaller 
than the surfaces' separation L, and showed its sinusoidal 
variation with the plates' lateral displacement [32 . In 
this situation good agreement was found with calcula- 
tions based on the PFA as the condition d <^ L was sat- 
isfied. More recently, the same group has measured the 
lateral force for asymmetric sinusoidal profiles accompa- 
nied by exact calculations [33l [34] using the method de- 
veloped in [17 . 

An interesting perspective for the lateral Casimir force 
lies in the realization of non-contact rack-and-pinion de- 
vices or ratchets [35H37] and of non-contact gears [38l [39] 
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which would allow for a force or torque transmission be- 
tween corrugated surfaces without bringing them in di- 
rect contact with each other. Typical profiles studied in 
this context are sinusoidal. 

In the present paper we study the Casimir energy be- 
tween two nanostructured surfaces as a function of their 
lateral displacement from which the magnitude of the 
lateral force may be easily deduced. We use the scat- 
tering approach j4QH42] to handle a variety of different 
periodic profiles. We present a detailed comparison of 
different profiles, such as circular, elliptical, triangular, 
trapezoidal and compare the Casimir energies they gen- 
erate to those of the commonly used sinusoidal and rect- 
angular profiles. Profiles with a base angle larger than for 
the rectangular profile generate stronger Casimir energies 
than profiles having a smaller base angle. While failing 
in general, we find the PFA to be a good approximation 
for the latter profiles displaced by half a corrugation pe- 
riod. Our results for circular and elliptical profiles can 
be applied to model the dispersive interaction between 
periodic arrays of nanowires or nanotubes [43H45] . 



II. SCATTERING FORMALISM 

For parallel plates, the Casimir force may be written in 
terms of scalar reflection coefficients, but for non-planar 
surfaces the specular reflection coefficients must be re- 
placed by general reflection operators that describe the 
non-specular diffraction by the surfaces [41^. In the case 
of two dielectric lamellar gratings, the Casimir interac- 
tion energy per unit area at zero temperature is given by 

nil 

E=J^ ///d^k^c[eindet(l-A^) (1) 

The Casimir energy is written in terms of reflection op- 
erators 1Zi and IZ2 which describe the diffraction by the 
two lamellar gratings. The operator e~'^^ accounts for a 
one-way propagation along the distance L separating the 
two gratings, with k, = ^^^^ + the imaginary longi- 
tudinal wave- vector and the transverse wave- vector, 
where all quantities are written at imaginary frequencies 
uj = after a Wick rotation. The operator M. thus 
represents one round-trip propagation between the two 
surfaces. 

Here we consider two dielectric lamellar gratings of ar- 
bitrary but symmetric profiles with a period d and corru- 
gation depth a separated by a vacuum slit such as shown 
in Figjl] In [17 a formalism has been developed to cal- 
culate the Casimir energy in this configuration for rect- 
angular profiles. We generalize this approach to arbi- 
trary profiles by dividing each corrugation line into K 
horizontal slices vertically stacked on each other. For 
a corrugation depth a, each slice is treated as a lamel- 
lar rectangular grating whose height along the axis is 




® z 



FIG. 1. Two nanostructured surfaces of arbitrary but sym- 
metric profiles with a period d and corrugation depth a are di- 
vided into K vertical stacks of rectangular slices (here K = 4). 

fixed dit a/K and whose length along the x— axis is given 
by the profile's length at the level of the slice. 

Let us first recall the main steps of the calculation for 
rectangular gratings separated by a vacuum slit. Because 
of time and z-invariance, we can write the electric and 
magnetic fields for i = (x, z) such as : 

Ei{x, y, z, t) = Ei{x, y) exp{ikzZ - iujt) (2) 

Hi{x, y, z, t) = Hi{x, y) ex.p{ikzZ - iujt). (3) 

In the following we will take /io = = 1, and c = 1 in 
vacuum. We need to find the longitudinal components 
outside the corrugated region {y > a) and within the 
transmitted region {y ^ 0). In the case of two planar 
interfaces, we can write the components of the fields 
in the vacuum region as an incident and reflected field: 

(x, y) = pe'^-^-'^yy + r^e^^-^+^^^^ (4) 

H,{X, y) = I^e^k^^-ikyy ^ ^h^ik^x+ikyy 

while in the material region, we can write the 
components of the fields as transmitted fields: 

E,{x,y)=t''e'^-''-'Ky (6) 

H,{x,y)=t^e'^-^-<y (7) 

where ky = u'^ — — k^ and ky = ecj^ — k^. — kl are 
the longitudinal wave-vectors in vacuum and inside the 
material, kx^kz are the components of the transverse 
wave-vector, e is the material's frequency dependent per- 
mittivity. The coefficients r^'^, t^'^ are respectively 
the incident amplitudes and the Fresnel- Stokes reflection 
and transmission amplitudes in the (e, /i)-basis of polar- 
izations, which are defined by imposing Hy = for the 
e-polarization and Ey = for the /i-polarization. 
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We may now generalize the above equations for grat- 
ings. Outside the corrugations {y > a), the fields are 
written as a Rayleigh expansion involving incident and 
refiected fields of order p and n respectively : 



Ez{x, y)y>a = Ip exp{iapX - ifS^^^y) 



(8) 



nez 



Hz{x, y)y>a = Ip ex.p{iapX - ip^^y) 

+ ^np exp(ianx + ip^^^y) (9) 



For the region inside the material ^ 0), the transmit- 
ted fields are given by : 

Ez{x, y)y^o = ^np exp(ianx - ip^^^y) (10) 

Hz{x,y)y^o = ^T^peyii[){ianX-il3^^y), (11) 
nez 

where we have used: 



ap 


— '^x ^ 


2iTp/d 


(12) 


an 
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2Tin/d 


(13) 




= u^- 
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(14) 
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(15) 
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(16) 



Ipi Rnp-, and Tnp are now the incidence, refiection, and 
transmission matrix elements respectively, n = cor- 
responds to a specular refiection. By symmetry, the 
other field components of the electric and magnetic fields 
can each be expressed through the z-components of both 
fields, following Maxwell's equations. 

We now need to determine the refiection coefficients 
Rnp of the rectangular corrugated gratings. To this aim 
we first rewrite Maxwell's equations inside the corrugated 
region {) < y < a through the set of first-order differential 
equations dyY = MF, for = {Ex, Ez, Hx, Hz) and 
M a constant square matrix of dimension 8A^ + 4. The 
solution for the fields is then of the form : 



F(^) = e^^F(O) 



(17) 



with : 



M 






ikzOln 



\ 
















(18) 



where the elements appearing in matrix M are block ma- 
trices of dimension 2N We can write the fields inside 
the corrugation region and match them through continu- 
ity relations for each Ex, E^, Hx, Hz, with equation (17), 



at boundary y = a ior y > a, and at boundary y = for 
y ^ 0. This allows to find the vectors F(a) and F(0), 
which can be written as a product of a matrix and the 
vector of variables X : 



F(a) = TX + y and F(0) = SX 



(19) 



with X^ = {R'^p,R^p,T^p,T^p,...), and Y being the 
variable-independent term including the polarization of 
the incident waves, since we must take into account the 



two polarizations, e and h independently: ±p 



(e) 



1 and 



Ah) 

Ip 

Ah) 

Ip 



for electric waves {Hz = 0), and 4"^^ = and 

1 for magnetic waves {Ez =0). Y hence char- 
acterizes the two separate solutions for e- and /i- waves. 
The solution is then of the form: 



X = {e'^'^'S - T)Y. 



We have 



x(/(^) = i,/^^) = o) = 



x(/^^) = o,/(^) = i 



/ p(e,e) 
' -Ltnp 
T-,(/i,e) 



/ -ttnp \ 
j^ih,h) 
rinp 



(20) 



(21) 



(22) 



so that we obtain the refiection matrix for each grating: 

-Llfip 



R{uj) 



iLnp 



-n{e,h) 
JXnp 



j^ih,h) 
iLnp 



(23) 



After making use of Cauchy's argument principle and 
normalizing the frequency by c, we arrive to the exact 
expression of the Casimir energy between the two grat- 
ings on a unit cell of period d and unit length in the 
^/-direction, with Ri{iC) and i?2(^0- 



E 



hdc 
8^ 



f^Tv/d 

d^ I dkz I dkx 



X lndet[l - Ri{ii)e-^^ R2{ii)e 



(24) 



with JC = diag(^^2 + A:^ + [kx + {2m7T/d)]^) and m = 

— X, . . . , -\-N. For the sake of clarity we have explicitly 
re-introduced the speed of light c here. 

III. ARBITRARY PROFILES 

We will now consider gratings of arbitrary symmetric 
profiles. The difference with the rectangular gratings ap- 
pears in the parameter di, which will now depend on y. 
Arbitrary profiles defined by di{y) can be divided into 
K slices, each of rectangular shape, as described in Fig. 



4 



(a) 



(triangular) 



^^ ^^^^ ^ 



(b) 



ttU 



(c) 



FIG. 2. Two-dimensional perspective on the different con- 
sidered periodic gratings, (a) and (c) are at the same scale, 
whereas the scale of (b) has been increased by a factor two. 



[T] For each slice (z), the spacing between the corrugation 
ridges is d^^ and the former scattering formalism for rect- 
angular corrugations can be applied. More specifically, a 
differential equation = M'^^^F akin to equation (17) 



can be solved within each slice {i) to relate the fields at 
boundary y = i-^ and y = (i + 1)-^. 

In a similar way than for the case K = 1 above, the 
field at y = a is thus related to the field ai y = via the 
relation : 



F(a) 



li=K 



F(0) 



(25) 



where the product Yl I'uns from i = K to i = 1. 

Hence a correct parametrization of the quantity di as 
a function of y allows one to generate arbitrary sym- 
metric profiles for the corrugations. The profiles that 
we will study as examples in the following are shown 
on Fig. |2]as cross sections while on Fig. |3]the grating 
structure becomes more apparent. Triangular profiles 
with di = d diS seen on Fig. [2k and Fig. |3] are gener- 
ated by the function di{y) =^^, sinusoidal profiles by 
di{y) = ((i/7r) arccos[l — (2?//a)]. We also study two types 
of trapezoidal profiles, having a base angle smaller than 
90° (Fig. |2|l) on the one hand, and having a base angle 
larger than 90° (Fig. ^p) on the other hand. They are 
characterized by djj < di and d^ > di respectively and are 
both generated by the function di{y) = {di — d}))y /a-\-d}). 
Ellipsoid profiles along the x- or ?/-axis are generated by 




(rectangular) 




(sinusoidal) 




(circular) 




(trapezoidal) 




FIG. 3. Three-dimensional perspective on the triangular, 
rectangular, sinusoidal, circular, and trapezoidal periodic pro- 
files shown in Fig. |2] The profiles are here represented at the 
same scale. 



di(^) ^ d-^-f^T^- {y-YY (Fig. and d^^y) = 
d — '^^R^ — {y — respectively, for R and r being 
the major and minor axes of the ellipse, and Y being the 
value of the ^/-coordinate of the ellipse center. Circular 
periodic profiles are also generated by these expressions 
(as seen on Fig. ^ and Fig. |3|, with R = r being the 
radius in the xy-plane. Note that what we call for ex- 
ample a circular or elliptical periodic profile is in fact a 
geometry where equally spaced parallel wires of circular 
or elliptical cross- sect ion cover the surfaces, as shown on 
Fig.[l 

Obviously a given profile will be better fitted for 
greater numbers of slices K. The number of slices thus 
determines the accuracy of the overall model. Fig. [4] 
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Number of slices K Lateral displacement (5/d 



FIG. 4. Dependence of the Casimir energy on the number 
of sUces K for two triangular gratings, with L = 100 nm, 
d 400 nm, a = 50 nm, and di{y) ^ + 200. 



shows the Casimir energy as a function of the num- 
ber of slices for two triangular profiles at a distance 
L — 100 nm, for a grating period d = 400 nm, corruga- 
tion depth a = 50 nm, and distance between the ridges 
di{y) = Ay -\- 200. For this example, convergence sets in 
for K ^ 20. 

In the following we will also compare the results of the 
scattering theory presented here with the proximity force 
approximation [14]. The PFA comes from the weighted 
sum of the planar normal contributions Epp{L) depend- 
ing on the local distances L within each period, and hence 
also on the lateral displacement 5 between the gratings. 
If we express the shapes of the arbitrary periodic gratings 
in an analytical form such as = f{x,6) for the lower 
grating and y = L-\-2a — f{x^5 = 0) for the upper grating 
in the xy-plane shown on Fig. [l] we can then define the 
function h{x, 6) = L+2a — /(x, S = 0) — /(x, S) in order to 
express the local distance of separation between the two 
profiles. Dividing the period d in a number N ^ oc of 
intervals of individual widths d/N 0, we then obtain a 
general expression of the Casimir energy in the PFA for 
arbitrary gratings as a function of lateral displacement S 



E^^^{L,S) = l r Epp{h{x,S))dx 
d Jo 

1 ^ / 



(26) 



x = z-.S 



IV. CASIMIR ENERGY FOR ARBITRARY 
PERIODIC GRATINGS 



We now evaluate numerically the Casimir energy for 
several types of profiles as a function of the surfaces' rel- 
ative lateral displacement S. The material chosen for 
these profiles is intrinsic silicon, which can be described 



FIG. 5. Dependence of the Casimir energy on the lateral 
displacement between two corrugated periodic profiles of in- 
trinsic silicon for L — 100 nm, K — 20, d — 400 nm, a = 50 
nm, and from top to bottom, di — Ay -\- 200 (triangular in 
black), = 3?/ + 200 (blue), di = 2y + 200 (red), sinusoidal 
profile (purple) ^ di — y-\- 200 (green) and di — 200 (rectangu- 
lar in orange). The sinusoidal profile is plotted for the same 
parameters, except di = (d/vr) arccos[l — {2y/a)\. 



by a Drude-Lorentz function 



(eo — eoo)^o 



(27) 



The numerical values are determined by realizing that 
at low-frequencies the dielectric function of intrinsic sil- 
icon approaches the constant value eo = 11.87 while 
with increasing frequency it is nearly constant and falls 
off only for high frequencies above a cut-off frequency 
ujQ ~ 6.6 • 10^^ rad/s. For high frequencies it reaches the 
asymptotic value e^o = 1.035 [47]. 

We first study the transition from triangular over 
trapezoidal to rectangular profiles, such as sketched by 
the colored shapes of Fig. [2^. Those profiles are 
parametrized from top to bottom by di {y) = Ay -\- 200 
(triangular: black), di{y) = 3?/ + 200 (trapezoidal: blue), 
di{y) = 2?/ + 200 (trapezoidal: red), di{y) = y ^ 200 
(trapezoidal: green), and di{y) = 200 (rectangular: or- 
ange) respectively. For these profiles the top distance 
di decreases successively from di = d = 400 nm to 
di = db = 200 nm by steps of 50 nm. The different 
Casimir energies for L = 100 nm, d = 400 nm, a = 50 
nm, and a number of slices K = 20^ are depicted in Fig. 
|5j Clearly the Casimir energy increases as di decreases, 
and this is especially true at 5 = d/2, meaning that both 
the Casimir energy and its modulation over lateral dis- 
placement are larger for smaller d/, i.e. for rectangular 
gratings and smaller for larger di such as triangular pro- 
files. Note also the behavior of the sinusoidal profile at 
S = (i/2, which shows that such profiles are much less 
sensitive to lateral displacement than triangular profiles 
or rectangular gratings. 

On Fig. |6] we show a comparison between the exact re- 
sults from the scattering formalism (solid lines) and the 
predictions from PFA (dashed lines) for the triangular. 
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0.0 0.1 0.2 0.3 0.4 0.5 

Lateral displacement S/d 



FIG. 6. Comparison between the exact results from scat- 
tering theory (solid lines) and PFA predictions (dashed lines) 
for the triangular (black) , sinusoidal (purple) , and rectangular 
(orange) profiles of Fig. [s] 

sinusoidal, and rectangular profiles (from top to bottom). 
Regardless of the considered profile, the PFA fails to cor- 
rectly describe the situation of no lateral displacement 
^ = 0. The error ratio £^pfa — ^scattering /^p fa there is 
approximately equal to 11% for triangular, 12% for sinu- 
soidal, and 15% for rectangular gratings. As one shifts 
6 to half-a-period S = the relative error for rect- 
angular gratings passes through zero to reach basically 
the same value of opposite sign at S = d/2. For rect- 
angular gratings PFA thus underestimates the Casimir 
energy when corrugation maxima face maxima and over- 
estimates it when maxima face minima. However, PFA 
turns out to give valid predictions for triangular and si- 
nusoidal gratings when they have a relative lateral shift 
of half a period. This is due to the fact that, unlike for 
rectangular gratings, in this situation the two triangular 
or sinusoidal profiles become parallel to each other. 

We now study the particular case of trapezoidal pro- 
files for which > di, with a base angle larger than 
rectangular gratings. We illustrate the results on Fig. [7| 
for L = 20 nm, d = 100 nm, a = 50 nm, and a num- 
ber of slices K = 20. The profiles are enlarged on the 
top while keeping the base constant via a parametriza- 
tion di = —O.by + 85 (grey), di = —O.by + 65 (deep 
purple), di = —0.5y + 45 (brown) from top to bottom. 
The dotted orange line corresponds to the Casimir en- 
ergy between rectangular gratings with di = di = 60 nm 
so that the corrugations are as wide as the top of the 
trapezoidal profile di = 85 — y/2. The difference in the 
Casimir energies associated with the rectangular profiles 
and the trapezoidal profiles, which have the same surface 
exposed in the near field, implies the existence of highly 
non-trivial mode contributions in the vicinity of the bases 
of the trapezoidal gratings. 

Finally we consider the Casimir interaction between 
two elliptical and two circular profiles, as shown in Fig. 
[Sj for 1/ = 100 nm, d = 400 nm, a = 50 nm, and a num- 
ber of slices = 15. The elliptical profile (lower curve 
in brown) has a major axis R = 50 nm along x, a minor 




-5000 H 



0.0 0.1 0.2 0.3 0.4 0.5 

Lateral displacement 6/d 

FIG. 7. Dependence of the Casimir energy on the rela- 
tive lateral displacement of two trapezoidal periodic profiles 
of intrinsic silicon for L = 20 nm, K = 20, d = 100 nm, 
a = 50 nm, and, from top to bottom, di = — 0.52/ + 85 (grey), 
di = —0.5^ + 65 (deep purple), di = — 0.52/ + 45 (brown). This 
is compared to rectangular profiles with di = di = 60 nm (or- 
ange dotted line). While we varied the spacing between the 
corrugations di, the grating period d was kept constant. 




0.0 0.1 0.2 0.3 0.4 0.5 



Lateral displacement S/d 

FIG. 8. Dependence of the Casimir energy on the relative 
lateral displacement between two periodic profiles shaped as 
circular (light blue) and ellipsoid (brown) from top to bottom 
for L = 100 nm, K — 15, d = 400 nm, a = 50 nm, and 
di 400 - 2y/50y - y'^ (circular with R = r = Y = 25 nm), 
and di — 400 — 4-^/50^/ — y'^ (ellipsoid with r = y = 25 nm 
and R — 50 nm). 



axis r = 25 nm, and its origin y = Y = 25 nm, so that 
di{y) = 400 — A^50y — y'^. The circular profile (upper 
curve in blue) has a radius R = 25 nm and its center 
at = r = 25 nm, such that di 400 - 2y/50y-y'^. 
The Casimir energy increases with the ratio R/r. This 
is again especially true at (5 = d/2. We may also note 
that compared to the profiles studied in Fig. |5) the en- 
ergy varies much more rapidly over S for the ellipsoid 
and circular profiles. This could be a consequence of the 
concave nature of these shapes for y < Y and could po- 
tentially affect the lateral Casimir force applications men- 
tioned earlier [35-39J. At fixed geometrical parameters 
and distances, the Casimir energy is overall much smaller 
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for these profiles than for the triangular and trapezoidal 
shapes discussed before. 

V. CONCLUSION 

We have studied the dependence of the Casimir energy 
on the lateral displacement for different arbitrary peri- 
odic gratings, ranging from triangular and sinusoidal pro- 
files, to trapezoidal, circular and ellipsoid shapes. Con- 
cerning the trapezoidal profiles, we find that at the same 
distance L, grating period (i, and corrugation depth a, the 
Casimir energy and its lateral modulation increase from 
triangular profiles over rectangular ones to those whose 
base angle is larger than 90°. This seems a consequence 
of the fact that for a given surface, the Casimir energy 
increases when the exposed surfaces in near-field is in- 
creased. Profiles with large base angles or at least rect- 



angular profiles seem therefore more promising for lateral 
Casimir force, Casimir torque or other non-contact de- 
vices |35H39] than sinusoidal or triangular profiles. How- 
ever, the exposed surface in near-field is not sufficient 
to estimate reliably the magnitude of the Casimir in- 
teraction even for small corrugation depth (a < L), as 
shown by the comparison of the Casimir energies between 
rectangular gratings and large base trapezoidal profiles. 
While failing to describe the Casimir interaction correctly 
in general, interestingly PFA gives valid predictions for 
sinusoidal and triangular gratings when they are rela- 
tively displaced by half a grating period. An interesting 
topic to further investigate would be to parametrize the 
profiles such that the lateral displacement 6 also depends 
on thus generating asymmetric profiles [34]. 

We thank the European Science Foundation (ESF) 
within the activity New Trends and Applications of the 
Casimir Effect (www.casimir-network.com) for support. 
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